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ABSTRACT 


The chromophores responsible for coloring the jovian atmosphere are embed- 
ded within Jupiter’s vertical aerosol structure. Sunlight propagates through 
this vertical distribution of aerosol particles, whose colors are defined by mo(A), 
and we remotely observe the culmination of the radiative transfer as I/F( A). 
In this study, we employed a radiative transfer code to retrieve w 0 {X) for 
particles in Jupiter’s tropospheric haze at seven wavelengths in the near-UV 
and visible regimes. The data consisted of images of the 2008 passage of Oval 
BA to the south of the Great Red Spot obtained by the Wide Field Plan- 
etary Camera 2 on-board the Hubble Space Telescope. We present derived 
particle colors for locations that were selected from 14 weather regions, which 
spanned a large range of observed colors. All mo (A) curves were absorbing 
in the blue, and mo(A) increased monotonically to approximately unity as 
wavelength increased. We found accurate fits to all m 0 (A) curves using an 
empirically derived functional form: cu 0 (A) = 1 — Aexp(— BX). The best-fit 
parameters for the mean mo (A) curve were A = 25.4 and B — 0.0149 for 
A in units of nm. We performed a principal component analysis (PCA) on 
our mo(A) results and found that one or two independent chromophores were 
sufficient to produce the variations in mo(A). A PCA of I/F( A) for the same 
jovian locations resulted in principal components (PCs) with roughly the same 
variances as the mo(A) PCA, but they did not result in a one-to-one mapping 
of PC amplitudes between the mo (A) PCA and I/F( A) PCA. We suggest that 
statistical analyses performed on I/F( A) image cubes have limited applica- 
bility to the characterization of chromophores in the jovian atmosphere due 
to the sensitivity of I/F{ A) to horizontal variations in the vertical aerosol 
distribution. 
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1 Background 


1.1 Chromophores and Two Ways to Define Color 

The jovian atmosphere displays significant contrast in its visible coloration. 
One or more coloring agent, or chromophore, is necessary to account for the 
observed color variations between weather regions such as the belts, zones, and 
vortices. However, the chemical identity, horizontal (latitudinal and longitudi- 
nal) and vertical distribution, and number of chromophores are still unknown 
(West et ah, 2004). 

When approaching the chromophore problem, we must make a distinction be- 
tween observed color and particle color. Observed color, I/F{ A), is the shape of 
a reflectance spectrum, which is a culmination of many wavelength-dependent 
radiative processes in the observed body. Particle color, U7 0 (A), is the shape 
of the single scattering albedo spectrum, which is a fundamental property of 
the individual particles. 

The relationship between observed colors and particle colors in the jovian at- 
mosphere is not necessarily a one-to-one function. Observed color is equivalent 
to particle color only in the case wherein isolated chromophores are observed 
directly, i.e., without any radiative contributions from intervening or adjacent 
material and without multiple scattering within the chromophore layer. These 
conditions clearly do not exist in planetary atmospheres: observed I / F{\) is 
the end result of radiative processes occurring within an atmosphere over dis- 
tances of tens to hundreds of kilometers that span large variations of pressure, 
temperature, aerosol abundance, and chemical composition. Thus, cuo(A) must 
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be derived through radiative transfer modeling of the vertical structure. 


1.2 Chromophore Characteristics via I/F(X) 

In order to spectrally characterize Jupiter’s atmosphere, most previous studies 
have analyzed the horizontal variations of // F( A) (e.g. Owen and Terrile, 1981; 
Beebe and Hockey, 1986; West et ah, 1986; Thompson, 1990; Dyudina et ah, 
2001; Simon-Miller et ah, 2001a). The main contribution of I/F( A) studies to 
the chromophore problem is an analysis of the directly observable quantities — 
I /F(X) as a function of latitude, longitude, and viewing geometry — for which 
any hypothesized distribution of chromophore candidates must be able to 
account. Although such studies describe real spectral and spatial distributions 
of observed colors, we cannot assume that the spectral and spatial distributions 
of the particle colors are equivalent. To some extent, the I / F( A) variations that 
we study are physically caused by horizontal variations in the atmosphere, such 
as cloud heights, optical depths, particle radii, etc. These factors may cause 
much larger variations in I /F{ A) than are caused by variations in zuo(X). 

For example, Thompson (1990) used a cluster analysis to find classifications of 
locations based on color and albedo. A cluster analysis creates classifications 
of spectra such that the spectra within a group are mutually similar and the 
spectral differences between groups are as large as possible. Each spectrum in 
the data set is assigned to a single group. Thompson (1990) analyzed 0.5° xO. 5° 
resolution maps from Voyager 2 acquired in three wideband filters centered at 
431, 564, and 599 nm and one narrowband filter centered at 621 nm, which 
spanned a weak CH4 absorption band at 619 nm. They found “20 distinct 
and five tentative units” for classification. While it is true that different com- 
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binations of vertical aerosol structure, viewing geometry, and chromophores 
resulted in the spectral differences defined by their classification scheme, the 
inclusion of a CH4 absorption band certainly increased their sensitivity to ver- 
tical aerosol structure and viewing geometry at the expense of sensitivity to 
chromophores. Therefore, their classification units are not directly applicable 
to the number or distribution of chromophores. 

Even if the method is an indirect one, an important product of I /F{ A) analyses 
is an estimate of the number of chromophores that are necessary to reproduce 
the observed variations in I/F( A). Strictly speaking, these estimates are nei- 
ther an upper nor lower limit to the actual number of unique chromophores 
that exist in the jovian atmosphere. Unless all other factors that are relevant 
to radiative transfer are ruled out, it is possible that only one chromophore is 
necessary. Alternatively, more than one chromophore may be present such that 
they always coexist spatially with the same mixing ratio and are, therefore, 
indistinguishable. 

A statistical method that has been used to find the number of independent 
variables contributing to a data set is principal component analysis (PC A) 
(for formalism, see Preisendorfer and Mobley, 1988; Murtagh and Heck, 1987). 
Unlike cluster analysis, PCA does not separate spectra into discrete groups. 
For data containing many spectra, PCA finds (1) a set of orthogonal spectral 
shapes 1 , which are called principal components (PCs), and (2) the amplitude 
of each PC that is necessary to reproduce each spectrum (typically the mean- 

1 We used the term spectral shape here (and hereafter, when referring to results from 
PCA) because PCs contain both positive and negative coefficients, which cannot be 
physically interpreted as actual spectra. In a given PC, wavelengths with coefficients 
of the same sign are correlated and those with opposite signs are anti-correlated. 
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subtracted spectrum) in the data by a linear superposition of all PCs. This is 
a factorization that can be represented as: 

V lX = (WH) iX = j2W ia Hax, (1) 

a=l 

where V is an m x r matrix that contains m spectra with r wavelengths, i 
denotes an individual pixel, A denotes an individual wavelength, and a denotes 
an individual PC. The PCs are chosen so that the first PC is in the “direction” 
(in r-dimensional wavelength space) containing the greatest amount of the 
spectral variance in the data. All variance in the direction of the first PC 
is then removed from the data. The second PC lies in the direction of the 
largest spectral variance that remains, which will necessarily be orthogonal to 
the direction of the first PC. This process continues until the number of PCs 
equals r, at which point the PCs constitute a complete set of basis vectors for 
the r-dimensional wavelength space. 

The advantage of PCA over cluster analysis in determining the number of 
chromophores is that one can estimate the number of unique spectral shapes 
that contribute to the total variance of the data. This is accomplished by 
estimating how many PCs describe real signal in the data and not merely noise. 
Each PC contains roughly the same amount of uncorrelated noise. Because 
each successive PC accounts for less of the total variance, it contains a smaller 
ratio of (correlated) signal to (uncorrelated) noise. PCs that have a higher 
order than the number of statistically independent physical processes that are 
affecting the shapes of the spectra will describe only noise and will all contain 
approximately equal amounts of the total variance. 

There is a noteworthy deficiency in all PCA results. We mentioned above that 
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the spectral shape of each successive PC is constrained to be orthogonal to all 
PCs that account for a larger amount of the variance. Thus, the spectral shapes 
of PCs 2 -r are not accurate representations of the true shapes of spectral 
variations in the data. 

Dyudina et al. (2001) used PCA to spectrally characterize Jupiter’s clouds 
based on Galileo Solid State Imager (SSI) and Near Infrared Mapping Spec- 
trometer (NIMS) data at 26 wavelengths. With their substantial wavelength 
coverage, inclusion of moderate and strong methane (CH 4 ) absorption bands, 
and inclusion of thermal emission at 5 /am, they found that 91% of the spec- 
tral variance was contained in the first three PCs. They found one PC that 
was associated with a violet-absorbing chromophore in a separate PCA of the 
Great Red Spot (GRS) and its immediate surroundings, but this was the only 
location in their data that yielded a PC with an anti-correlation between vi- 
olet reflectivity and all other non-thermal bands. However, only one of their 
wavelengths (410 nm) was in the spectral range of the broad chromophore 
absorption blueward of ~600 nm (see Fig. 1), and their PCA results were 
mostly indicative of spatial variations in vertical aerosol structure, which is 
consistent with their interpretations of them. 

Simon-Miller et al. (2001a) conducted a PCA study that was focused on chro- 
mophore absorption. Their data consisted of I /F-calibrated global maps of 
Jupiter from HST images acquired in October 1995 and October 1996. The 
first set contained continuum filters F410M, F555W, and F953N, and the sec- 
ond set added the F673N filter. They determined that only three spectral 
components are required to explain the deviations from the mean I jF{\). 
PCI, containing 91% of the variance, was spectrally gray, and therefore did 
not correspond to a chromophore. They concluded that PCI “is probably re- 
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lated to density but is not exactly cloud optical depth”. PC2 contained 8% of 
the variance and indicated the presence of a red chromophore. PC3, presum- 
ably arising from the presence of a second chromophore, contained 1% of the 
variance and was present in the GRS and some smaller anticyclonic ovals. The 
1996 data set with four filters yielded a fourth PC with a variance <C 0.1%, 
which was less than the observational error. 

Thus, Simon-Miller et al. (2001a) constrained the number of statistically in- 
dependent factors affecting their filter set to three, attributed two of the three 
to chromophores, obtained a map of each PC’s horizontal distribution, and 
provided a low-resolution spectral shape — though affected by the orthogo- 
nality constraint — for the most prominent I / F{ A) color variation due to a 
chromophore (PC2). But to conclusively identify a chromophore, we need to 
do more than isolate differences in observed color through analyses of I / F{\). 
We must eventually have a spectrum of the particle color itself. 


1.3 Chromophore Characteristics via zuq(X) 


Retrievals of vjq(\) must be made in the context of an atmospheric model. The 
complications and degeneracies inherent in radiative transfer models explain 
why fewer studies have set out to characterize rt7o(A) than I/F( A). The goal 
of many radiative transfer studies is to explore the vertical distribution of 
aerosols, and their data do not contain the appropriate spectral coverage to 
study rao(A). It is often impractical to obtain center-to-limb feature tracks 
in several — preferably narrowband — continuum filters in addition to the CH4 
absorption band filters that are required for vertical discrimination. 
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Simon-Miller et al. (2001b) studied the vertical aerosol structure and particle 
absorption in the continuum with a radiative transfer analysis of Galileo SSI 
data. 410 nm and 756 nm were the only continuum wavelength for which 
they had data. They fixed ro 0 (756)=l and explored the range of cuo(410) — 
not the number of chromophores nor their spectra. They preformed retrievals 
of vertical aerosol structure and cu o (410) for zones, belts, vortices, and a hot 
spot. As in the Banfield et al. (1998b) study, which Simon-Miller et al. (2001b) 
extended by considering particle color, the vertical structure that Simon-Miller 
et al. (2001b) employed contained as few layers as possible to fit their data. For 
the majority of locations, their data were “well fit with a three-layer model 
that contains a thin stratospheric haze, a denser tropospheric haze, and a 
tropospheric cloud sheet at pressures consistent with an ammonia cloud.” 
The vertical distribution of their cuo(410) retrievals showed the chromophores 
to reside mostly in the tropospheric haze. 

Perez-Hoyos et al. (2009) presented a detailed radiative transfer study of 
Oval BA and its immediate surroundings to investigate its change in color 
from white to red (Simon-Miller et al., 2006). They conducted a radiative 
transfer analysis of HST images acquired before (2005) and after (2006) the 
color change, retrieving n7 0 (A) in six filters: F250W, F330W, F435W, F550M, 
F658N, and F892N. They also used a three layer model and found the chro- 
mophores to primarily reside in the tropospheric haze, and we adopted this 
hypothesis for the chromophore location in this work. Their zuq(\) retrievals 
for the Oval BA region are the most complete spectra of jovian chromophores 
of which we are aware. However, they are not necessarily descriptive of the 
global chromophore population because it was retrieved from an isolated re- 
gion of the planet. 
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Our goal in this work is to use sufficient spatial and spectral coverage to char- 
acterize the particle colors of the global chromophore population. We retrieved 
ruo(A) from multispectral HST WFPC2 images and statistically analyzed the 
horizontal variations. Section 2 describes the observations, data reduction, and 
data selection for modeling. Section 3 presents the radiative transfer code, the 
modeling procedure, and an estimate of the expected degeneracy. Section 4 
presents the retrievals and an analysis of the particle color variations. Section 5 
discusses the derived chromophore characteristics and explores the extent to 
which an I/F( A) analysis reaches similar results. Section 6 contains the con- 
clusions. 

2 Observations 

2. 1 Overview 

HST WFPC2 observations of Jupiter were made on 15 May, 28 June, and 8 
July 2008, covering a passage of Oval BA south of the Great Red Spot (GRS). 
Two separate epochs of data exist for 15 May due to the separation between 
transit times for the GRS and Oval BA. Nine filters sampled the continuum 
(Fig. 1), and seven of these were narrowband filters, making this data set ideal 
for color studies at high spatial resolution. 

We created a spectral image cube for each epoch of observation. These cubes 
contain two spatial dimensions (planetographic latitude and System III lon- 
gitude) and one spectral dimension. Associated with each mapped image in 
each cube are maps of the three relevant viewing geometry quantities: p, po, 
and 0, where po is the cosine of the incidence angle, p is the cosine of the 
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emission angle, and 0 is the phase angle. 

Each spectral image cube contains at least one image per filter. In the cases 
where two images exist for one filter, both images with their unique viewing 
geometry data were used. Two F410M images exist for both epochs on 15 
May, and two F673N images exist in all epochs. Table 1 contains a summary 
of the observations and the spectral image cube assignments. 


2.2 Absolute Photometric Calibration, Navigation, and Mapping 

We converted images to I/F using the PHOTFLAM photometry keyword 
provided by the Space Telescope Science Institute for each filter, the solar flux 
across each pass band, distance of the Sun to Jupiter, image integration time 
and pixel solid angle, as in Simon-Miller and Gierasch (2010). In addition to 
the standard conversion, we corrected the F255W filter throughput for long- 
term degradation following Gonzaga et al. (2006). After conversion to I/F, a 
comparison of the approximately full disk I /F(X) to a standard Jupiter full- 
disk spectrum (Karkoschka, 1998) showed that the F343N filter brightness 
values appeared anomalously low. This was subsequently confirmed via stan- 
dard star monitoring, and a new filter throughput was calculated, changing 
the final I/F by nearly a factor of two (see Gonzaga and Biretta, 2009). With 
this new value, the final full-disk I/F( A) is consistent with Karkoschka (1998) 
(Fig. 1). 

To properly determine planetary coordinates for mapping, we navigated each 
image. The navigation was performed using the known sub-spacecraft and 
sub-solar points, the camera’s plate scale and focal length, and the distance to 


13 


Jupiter to size an ellipse with the correct aspect ratio and terminator and limb 
orientations. The north pole angle was based on Jupiter’s position angle on the 
sky and approximate roll angle of the camera. An initial location was chosen 
for the ellipse center by manually matching the limb and terminator locations, 
usually to within a few pixels. An automated routine then located actual limb 
and terminator points by contrast with the sky background, searching within 
20 pixels of an initially estimated ellipse location. The planet center and north 
angle were then iterated using a least squares minimization until residuals were 
less than 0.1 pixels, and usually less than 0.05 pixels for most wavelengths. 

However, the precision of the navigation was better than the accuracy, which 
was revealed after we reprojected the images onto cylindrical maps using 
MaRC 2 . This software package enabled us to perform reprojections of I / F, fj,, /io , 
and (f) at each pixel in the HST images based on Jupiter’s oblate spheroidal 
shape and viewing geometry. With the cylindrical maps we created spectral 
image cubes for each epoch (Fig. 2) and visually inspected the relative spatial 
positions of all features across all maps in each cube. Although we attempted 
to solve minor registration issues using new planet navigations, this procedure 
did not yield satisfactory results. We found it necessary to manually apply 
a small shift (< 0.29° planetographic latitude, < 1.43° longitude) to indi- 
vidual maps in order to have the same latitude and longitude position for a 
given atmospheric feature in all maps within a given spectral image cube. The 
maximum displacement of features due to zonal flow during the time elapsed 
between the first and last images in each image cube was on the order of one 

2 MaRC (Map Reprojections and Conversions) is a free C++ library and program 
developed by Ossama Othman to produce map projections. Further information 
and software can be found at http://sourceforge.net/projects/marc/ . 
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pixel. We therefore did not attempt to correct for this effect. 


2.3 Data Selection 

Our goal in this work was to demonstrate the possible range of jovian chro- 
mophores’ particle colors. To achieve this, we needed to retrieve vjq(X) from 
jovian locations with significantly different I/F( A). As the data did not con- 
tain feature tracks from center to limb at any wavelength, we could not retrieve 
U 7 ()(A) for individual pixels. Rather, it was necessary to identify groups of pixels 
to model simultaneously (hereafter, data cuts) in order to improve the range 
of fi and /< 0 — thereby improving our confidence in the model retrievals — as 
much as possible. Each data cut is extracted from a set of locations that we 
expect to have similar vertical structures based on spectral similarity. It is not 
possible to construct a data cut for which it is known a priori that a single 
vertical structure is contributing to I/F( A). However, we do know a priori 
that a given vertical structure will produce identical I/F( A) at two locations 
as 5p,,5fi o, and S<p approach zero. The I / F(X) of a given vertical structure 
should have smooth center- to-limb variations because in this case the I /F( A) 
will only change as a function of the viewing geometry. Therefore, locations 
with vertical structures that are different from that of a fiducial location can 
be excluded (to first order) from a data cut by selecting only spectrally similar 
locations. 

To create a data cut, we first selected a fiducial pixel on the mapped data 
cube and included all pixels within ±0.39° latitude (±5 pixels in the mapped 
images) with similar spectra (for an example, see Fig. 3). The criterion for the 
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spectral similarity of a test pixel with the fiducial pixel was: 

|/ / P'o,fiducial{ A) // Fq t test (A)l nos 

I/F 0Jlducial (X) ' ’ 


for all values of A, where: 


I/F 0 = 


IJF' 

Mo 


(2) 


(3) 


Thus, the reflectivity of all pixels in a data cut were within 5% of the fiducial 
pixel’s reflectivity. The selected pixels were split into bins 1° wide in longitude. 
We averaged the I/F, /t, /jq , and <p for the pixels in each bin in each image of 
the data cube. 

The quantity //To was used to compare pixels because this eliminated spectral 
differences due solely to variable illumination. Doing so assumed a value of 
k — 1 in the Minnaert, function for limb darkening: 

A = A M ^ k ~ l , (4) 


where A is albedo and Am is the Minnaert albedo. However, we cannot as- 
certain the true value of k for a given location because of the lack of multiple 
images of the same feature as its viewing geometry changed. If our assumption 
that k = 1 is significantly in error, then there are two consequences for the 
data selection as the \i and //, 0 distances from the fiducial point increase: (1) 
pixels with identical vertical structure to the fiducial point may be wrongly 
excluded in the data cut because their I/Fq(X) now differs by >5%, and (2) 
pixels with different vertical structures from the fiducial point may be wrongly 
included in the data cut. 

Based on visual inspection of the data cuts produced by this method, it ap- 
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peared to satisfactorily select a uniform set locations similar to the fiducial 
point and reject locations with obviously different clouds. In this sense, it is 
an improvement upon the practice of including (in one data cut) all longitudes 
within a narrow latitude range despite the visible heterogeneity that exists in 
the jovian atmosphere. Our method did, however, appear to prevent our data 
cuts from including longitudes as close to the limb of the disk as would have 
otherwise been possible. But, this is of secondary importance compared with 
successfully excluding pixels that do not belong to the same spectral class as 
the fiducial point, which we believe that our method reasonably attains. 

We selected the fiducial pixels for data cuts from 14 jovian weather regions 
(Table 2). We did not perform any spectral analysis to determine which fiducial 
points to select. We visually identified zonal bands with smooth brightness 
variations in longitude and selected fiducial points manually from them. We 
did not include highly inhomogeneous bands, such as the South Equatorial 
Belt, or bands dominated by vortices, such as the southern zones. We also 
visually identified the vortices and discretely colored regions within the GRS 
and BA to manually select fiducial points. The spatial coverage of our data 
cuts is globally representative in the sense that zones, belts, and vortices of 
various colors are included. 

Figure 2 displays all locations used in this study marked on an approximate 
true-color composite of the data, and Figure 4 shows the average reflectance 
for all data cuts. Neither the modeling nor subsequent statistical analysis 
took the regional groupings into account — they are used only as a reference 
for ease and clarity of discussion. Six zonal regions were selected based on 
their visible uniformity in longitude, which indicated that many pixels could 
be grouped into a data cut appropriately: the Equatorial Zone (EZ), the North 
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Equatorial Belt (NEB), the North Tropical Zone (NTrZ), the southern area of 
the North Temperate Belt (NTBs), the northern area of the North Temperate 
Belt (NTBn), and the North Temperate Zone (NTZ). We also selected eight 
smaller-scale systems based on their spectral uniqueness: the red annulus, 
bright annulus, and red center of the GRS; the red annulus and bright center 
of Oval BA; a small red anticyclone that passed through the perimeter of the 
GRS (Oval 2); white oval A5; and a red oval embedded in the NEB (hereafter, 
NEB1). Figure 5 shows the average reflectance for all data cuts within each 
weather region studied. In total, we selected 1345 data cuts from these 14 
regions. 

As seen in Table 2, we used many more data cuts from the smaller regions 
than the zonal regions. The reasons for this are derived from the fact that the 
zonal regions contain many more samples per data cut. The greater the range 
of p and po that a data cut contains, the better the model fits are constrained. 
For longitudinally restricted weather systems such as the GRS, the data cuts 
have only a few p and po points, and these are closely spaced. For zones and 
belts, however, we had extensive coverage in p and po by including a large 
range of longitudes. In order to increase our confidence in model retrievals of 
the small systems, we modeled many more data cuts taken from them. 


3 Radiative Transfer Model 

3. 1 Model Heritage and Methodology 

The radiative transfer code we employed was originally developed by Banfield 
et al. (1998b) to analyze Galileo SSI data. Their adding and doubling code was 


18 


based on the work of Hansen and Travis (1974). Banfield et al. (1998b) utilized 
data taken in filters centered at the strong CH 4 band at 889 nm, the moderate 
CH 4 band at 727 nm, and the continuum at 756 nm. Images in each filter were 
selected for three viewing angles in which the SSI target location was near the 
terminator, nadir, and limb. The range of pressure levels probed by the three 
filters combined with observations from multiple viewing angles and small- 
scale spatial variations of I/F{ A) within each data cut enabled Banfield et al. 
(1998b) to retrieve optical depths and vertical positions of aerosol layers. 

Simon-Miller et al. (2001b) supplemented the code by including images from 
the SSI 410-nm filter and solving for wq at 410 nm in each aerosol layer. Our 
work extends the Simon-Miller et al. (2001b) model for use with HST WFPC2 
filters and an arbitrary number of images in each filter. For our study, we used 
Rayleigh gas that contained H 2 , He, and CH 4 with number mixing ratios 
with respect to H 2 of 0.157 for He (von Zahn et al., 1998) and 0.00237 for 
CH 4 (Wong et al., 2004). Following Banfield et al. (1998b), we assumed that 
the aerosols are spherical particles, thereby using Mie scattering theory to 
determine extinction efficiencies and phase functions. We assumed the aerosol 
composition to be ammonia ice with an index of refraction of 1.4 (Martonchik 
et al., 1984). We employed a gamma distribution of particle radii (Mishchenko 
et al., 1999): 

n(r) = c x r^~ 3b ^ b exp 

where n is the number of particles with radius r, ro is the typical particle ra- 
dius, b is a characteristic width set to 0.1, and c is set such that f!. max n(r) = 1. 
For the particle phase function, we used a Henyey-Greenstein phase function 
with a single asymmetry parameter fit to the Mie phase function for a given 


r 

r 0 b 


be (0,0.5), 


(5) 
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particle size distribution. We parametrized each aerosol layer by its base pres- 
sure, optical depth, and typical particle radius. 

There are three aerosol layers and twelve model parameters per layer in our 
code: base pressure (P), optical depth (r) at 502 nm, particle radius (r 0 ), 
and single scattering albedo (cu 0 ) for each of nine filters. Table 3 contains gas 
absorption parameters and Rayleigh scattering parameters. We calculated gas 
absorption by averaging the CfR absorption coefficients, k, from Karkoschka 
(1998) over the total system response curve of each filter. We fit both gas 
absorption and Rayleigh scattering with a single-parameter Beer’s Law. 

We began this work before Karkoschka and Tomasko (2010) published an up- 
dated model of CH4 &- values. However, the differences in the &- values between 
the two models were negligible for our filters. The &- values reported are so low 
that Rayleigh scattering entirely dominates gas absorption in the radiative 
transfer, and we continued with the k - values of Karkoschka (1998) due to its 
simpler treatment. The relative importance of CH4 absorption and Rayleigh 
scattering is illustrated in Fig. 6, which shows the two-way transmissivity at 
each wavelength due to CH4 absorption and Rayleigh scattering in a clear, 
aerosol-free atmosphere. Two-way transmissivity is the fraction of light that 
remains in the beam — i.e. the fraction that is neither absorbed nor scattered 
out of the line of sight — for a path from space down to a given pressure and 
back. The two-way transmissivity curves indicate the vertical sensitivity of the 
filters. The vertical sensitivity for this set of filters due to CH4 absorption is 
limited to pressures >10 bar. 

We estimated the uncertainty in relative photometry (a) in the ///-’-calibrated 
images by comparing the standard deviation of pixel values within each bin 
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across the weather regions. This is used to calculate the value of x 2 f° r each 
model: 




1 


g (» - /,) 2 


W'obs Tlpar & 


(6) 


where n 0 b s is the number of observations, n par is the number of free parameters, 
Hi is the observed I/F at the viewing geometry and wavelength of observation 
i, fi is the corresponding model result, and a is the standard deviation of the 
filter used for observation i. 


3.2 Vertical Structure and Parameter Space 

Figure 7 shows the format for the vertical structure used in our models. We 
included three adjoining aerosol layers to model pressures <1 bar, similar 
to the standard model presented in West et al. (2004) (see their Fig. 5.15): 
a stratospheric haze, tropospheric haze, and tropospheric cloud. We use the 
term “haze” to refer to aerosol layers that contain both Rayleigh-scattering 
gas and aerosols. We call a layer that solely contains aerosols a “cloud” . 

The base of the stratospheric haze is attached to the top of the tropospheric 
haze, and the base of the tropospheric haze is attached to a tropospheric cloud, 
which is located in the vicinity of the ammonia condensation level. The top 
and bottom of the cloud are set to the same pressure. The base of the model 
atmosphere is a semi-infinite Rayleigh-scattering gas layer. 

The infinitely-thin tropospheric cloud — which was also used by Banfield et al. 
(1998b) and Simon-Miller et al. (2001b) — is a numerical convenience and is not 
necessary to fit the data. It results in an aerosol layer without any Rayleigh- 
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scattering gas mixed in, which is a convenient approximation when aerosols 
dominate the radiative transfer. This approximation becomes increasingly de- 
generate in pressure as the total r of the true cloud decreases and as the 
vertical extent of the true cloud increases. For an optically thick cloud, the 
pressure level of the infinitely thin cloud is where the top of the physically 
extended cloud would be located. For an optically thin cloud, the pressure 
level of the infinitely thin cloud would probably be near the middle of the 
physically extended cloud. 

We implemented several parameter constraints to avoid fitting degenerate pa- 
rameter combinations. Table 4 shows all fixed and free model parameters, 
and Table 5 contains the initial parameter values. Due to the limited vertical 
discrimination possible with our data, we fixed the pressure parameters. To 
explore the relevant pressure range of tropospheric clouds, we used three dif- 
ferent models with tropospheric cloud pressures at 0.4, 0.7, and 1.0 bar. The 
interface between the stratospheric haze and tropospheric haze was fixed at 
the approximate location of the tropopause (0.1 bar) in all models. Although 
some studies have found that a clearing in the aerosols near the tropopause 
was necessary to fit their data (Banfield et ah, 1998a; de Pater et ah, 2010a), 
we do not believe that we can discern a difference between a gap or no gap in 
the aerosols with our data due to degeneracies with the pressure parameters 
(see §3.4). 

We fixed the particle radii at the central values of the ranges reported by 
Simon-Miller et al. (2001b) after finding very little ability to constrain radii 
retrievals with our data. We derived r for all three layers but constrained 
them such that r < 50. Without a clear ability to discern the altitude of each 
aerosol layer, we are not able to speculate in this work about the aspects of 
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the vertical structure or changes therein that are associated with chromophore 
production/ visibility. See Perez-Hoyos et al. (2009) and de Pater et al. (2010b) 
for discussions of chromophore variability mechanisms. 

We fit zuq (255-343 nm) only in the stratospheric haze due to the shallow pen- 
etration depth of UV light, and fit m 0 (375-673 nm) only in the tropospheric 
haze. Retrieved U7 0 (255) and tuo(343) values are shown in a few instances but 
are not discussed. They are necessary parameters for the radiative transfer 
model, but they are not well constrained. 

We tied the wq (255-343 nm) in the tropospheric haze to the value of cuo(375 
nm) for several reasons. The sensitivity of the model to u 7 0 ( 255-343 nm) in the 
tropospheric haze was very low. When we left these as free parameters while 
testing the model, they often converged to values near zero. This increased 
the values derived for cuo (255-343 nm) in the stratospheric haze, which must 
be free parameters to properly fit T\ in the stratosphere. We concluded that a 
physically reasonable solution was to force cu 0 (255-375 nm) in the tropospheric 
haze to have a spectral slope of zero. We assumed conservative scattering 
(w o=1.0) for all fixed z&o parameters. 

Simon-Miller et al. (2001b) determined that Jupiter’s I/F at 410 nm was 
consistent with models in which chromophores were located only in the tropo- 
spheric haze above the main ammonia cloud deck. We restricted the aerosol 
particles’ visible coloration to the tropospheric haze based on these results. 
Although Smith and Tomasko (1984) found the coloration to reside in both 
the tropospheric haze and the main ammonia cloud deck, we limited our pa- 
rameter space here in an effort to find the simplest model that can sufficiently 
reproduce the data. 
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3.3 Parameter Space Search 


The original y 1 2 convergence algorithm used by Simon-Miller et al. (2001b) 
can fit two to four parameters at a time. We found that more than four free 
variables at once will undermine the fidelity of the minimization. Simon-Miller 
et al. (2001b) had twelve parameters to fit (four per model layer), and they 
performed their parameter space search manually. The user selected a small 
number of free parameters and ran the y 2 convergence algorithm. The user 
then iterated this process, choosing a new set of free parameters each time, 
until the user judged that the model provided a satisfactory fit to the data. The 
parameter space in this work also included twelve free parameters, however 
the number of data cuts to be modeled was large enough to warrant using an 
automated search algorithm in conjunction with the y 2 convergence algorithm. 
This operated in three phases. Phase A determined the appropriate range 
of values for each parameter. This defined the working search space. Phase 
B conducted a grid search of this space, and Phase C converged on a final 
solution. 

For Phase A, we assumed an initial vertical structure model (see §3.2) and fit 
two free parameters at a time by an iterative process: 

(1) Estimate the increment in the model parameters that would minimize y 2 
(Eq. 6). We used a singular value decomposition on a design matrix con- 
sisting of the difference between the observations and the current model 
predictions along with the partial derivatives of the observables with re- 
spect to each model parameter (Press, 2002). 

(2) Change the model parameters by 90% of the increments estimated in 
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step 1. (Moving 100% of the distance may overshoot the actual minimum, 
but choosing to move too slowly would be computationally expensive.) 

(3) Run the radiative transfer code on the new model parameters and eval- 
uate y 2 . 

(4) Repeat steps 1 and 2 until convergence is reached, and save the final 
parameter values. 

(5) Reset the model to the initial values, and proceed to the next parameter 
pair. 

Each parameter was paired more than once, and the value to which a pa- 
rameter converged depended upon which other parameter was free to vary. 
For example, wq(673 nm) in the tropospheric haze may converge to one value 
when varying simultaneously with the optical depth of the stratospheric haze 
and a different value when varying with the optical depth of the tropospheric 
haze. The maximum and minimum values to which each parameter converged 
in Phase A defined the parameter space for the grid search in Phase B. We 
explored every free parameter with at least two pairings, but used just 20% 
of all possible pairings. Using them all would have greatly increased the re- 
quired computing time. Run times for each phase on a Pentium 4 processor @ 
3.0 GHz ranged from a few minutes to a few hours, depending on the number 
of samples in the data cut. 

In Phase B, we generated model results for a parameter grid that covered the 
Phase A range. If the range of a given parameter was <1%, we then used 
the average value only. The maximum number of grid points, n g , was set to 
n g, m ax = 9 99 in order to prevent grids from becoming unmanageably large. If 
Phase A returned n g > n g _ max , then the grid resolution for one parameter was 
reduced. This process was iterated until n g < n gmax . The y 2 was evaluated 
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for all models in the resulting grid. If the best-fit model did not meet the y 2 
cutoff (empirically set to 2.0), we then set this model as the initial conditions 
for another iteration of Phases A and B. 

If the minimum y 2 of the grid was < 2.0, we then proceeded to Phase C. 
In this last phase, we ran the parameter-pairing method described in Phase 
A but without resetting the parameters to the initial values between each 
iteration. This resulted in a final model in which the parameter values were 
not quantized, unlike the parameter values in the grid from Phase B. 


3-4 Tests of Methodology and Parameter Assumptions 

An example of the best-fit model to a data cut from the EZ using our method- 
ology (Fig. 8) shows that we find appropriate fits to I/F considering both ge- 
ometry and wavelength. The parameter space explored by the y 2 - minimization 
algorithm in Phases A and B and the best-fit model are shown in Figure 9. 
It can be seen that the values checked for each parameter extend above and 
below the best-fit value. 

We wanted to know whether or not our parameter search finds an appropri- 
ate y 2 minimum using the initial parameters that we selected. We therefore 
created a set of synthetic data cuts using parameters that span most of the 
parameter space of the retrievals from the real data presented in Section 4. 
We used one data cut per epoch from each of the six zonal regions (EZ, NEB, 
NTrZ, NTBs, NTB 11 , NTZ), which is a total of 24 data cuts. We took the 
best-fit model result (i.e. synthetic data generated by the radiative transfer 
code) from each of the three models (P 3 = [0.4, 0.7, 1.0] bar) and added Gaus- 
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sian noise (a = 0.005 [/ / F] ) to visually match the amount of scatter present 
in the real data. We retrieved parameters for these synthetic data cuts with 
the model with the corresponding value of P3. For example, our modeling of 
the real data cut from the EZ in Cube A resulted in three best-fit synthetic 
data cuts, one for each tropospheric cloud pressure, P 3 . We then added noise 
to each of these results. We fit a model with P 3 = 0.4 bar for the synthetic 
result with P3 = 0.4 bar, a model with P 3 = 0.7 bar on the synthetic result 
for P 3 = 0.7 bar, and a model with P 3 = 1.0 bar on the synthetic result 
for P3 = 1.0 bar. The retrievals (Fig. 10) demonstrate the ability of the y 2 - 
minimization algorithm to find the correct local minimum (though it may not 
be a truly global minimum) in the parameter phase space when starting with 
the initial parameters used for the data retrievals. 

One of our simplifying assumptions was to leave the pressure of the strato- 
spheric/tropospheric haze interface (Pi) fixed at 0.1 bar. In reality, this inter- 
face undoubtedly varies in height across the planet. Banfield et al. (1998b), 
Simon-Miller et al. (2001b), and Perez-Hoyos et al. (2009) report retrievals 
in the approximate range of 0. 1-0.3 bar. To estimate the effect on se7 0 (A) re- 
trievals due to fixing this parameter value, we generated synthetic data with 
Pi=0.2 bar and performed retrievals under our assumption of = 0 . 1 bar. 

We created this synthetic data differently from that used in the previous test. 
We selected the data cut with the most samples from each weather region, 
changed the best-fit model with P3 — 0.7 bar so that P\ = 0.2 bar, and ran 
the radiative code on this altered model to generate the synthetic data. In this 
test, we used all three models to fit the data generated with P 3 = 0.7 bar. 

Figure 11 shows the retrieved parameters versus the true parameters, and 
Figure 12 shows the true tuo(A) and retrieved zvo(A) for the tropospheric haze 
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particles from all three models by weather region. Prom Figure 11 we can 
see that the retrieved zu 0 (X) values differ at most by ~10% at 375 and 390 
nm but the differences quickly decrease with longer wavelengths. The largest 
wo(X) differences are seen in the synthetic data with the darkest particles. 
The optical depths retrieved differ by as much as about a factor of five for t\ 
and about a factor of two for r 2 and r 3 . Also note in Figure 12 that in general 
the three retrieved models (gray lines) have significantly less variance between 
themselves than the difference between the true tu 0 (X) (black line) and the 
retrieved ones. This demonstrates that the dependence of the z&o(X) retrievals 
on P 3 is negligible compared to Pi in all regions except the EZ. 

Another set of simplifying assumptions was to fix the typical particle radii in 
each of the three aerosol layers. The particle radii in our models determine 
the scattering efficiency of aerosols, which is normalized by the scattering 
efficiency at 502 nm. The optical depths that we report are at this wavelength. 
If the actual particle radii differ from our assumed distributions, then our 
assumed r( X) relationship will be incorrect. One can compensate for the over- 
or underestimation of the model reflectivity due to an incorrect r( A) curve 
by adjusting mo(A) and still produce satisfactory fits to the data. This is a 
significant case of degeneracy in radiative transfer models. 

To explore this degeneracy, we created synthetic data tests in the same man- 
ner as the Pi test above. We conducted three separate tests to separate the 
effects of changing the particle distribution in each aerosol layer. We used our 
initial models (7*1 =0.03, r2=0.9, r 3 = 2.0) to fit data generated with ri=0.1 /an 
(Figs. 13 and 14), t’2=2.0 /urn (Figs. 15 and 16), and r 3 =4.0 /mi (Figs. 17 and 
18). Fitting with an incorrect particle size distribution had a similar effect 
for all three aerosol layers. The optical depth results for all three suffered a 


28 


similar level of inaccuracy as the Pi test. Also, as with fitting with an incor- 
rect Pi, the largest u7o(A) differences are seen in the synthetic data with the 
darkest particles and the blue end is worse than the red. The ri=0.1 yum test 
yielded the most disparate results, with the retrieved zu 0 (375) up to ~50% too 
low. However, the retrievals for g 7 O (410) were generally accurate to <6% and 
w 0 (437-673) to <3%. 


4 Results 


We derived best-fit models of r and u7o(A) for all 1345 data cuts. The data were 
well fit; the models typically had (reduced) y 2 values of ~1. Note that degen- 
eracy exists in all radiative transfer solutions to tropospheric structure (West 
et ah, 2004). Therefore, these results are indicative of the true structures — 
seen through the constraints of our simplified model — and are not to be taken 
as precise representations. 

The results we report for each data cut are the average of all forward models 
computed during the y 2 -minimization routine that have y 2 < 2yf nin , where 
Xmin IS the value for the best fit. Most data cuts were well fit by the parameter 
space search using all three initial models, which varied only in the pressure 
level of the interface between the tropospheric haze and cloud. Therefore, the 
results presented typically contain contributions from at least one model from 
each vertical structure. 
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4-1 Aerosol Optical Depths 


The particle color retrievals are dependent upon the optical depth retrievals. 
Much data exist concerning Jupiter’s vertical aerosol structure, and little data 
exist for the particle color — none of which describes the global distribution. 
Therefore, we first considered whether or not the vertical aerosol structure 
retrievals provided a realistic physical context for the particle color retrievals. 

The optical depths retrieved for each aerosol layer are presented in Table 6 
and Fig. 19-21 for each weather region. The number of data cuts per weather 
region is noted in each panel. The stratospheric haze (Fig. 19) has T\ < 0.35 
in all weather regions and a typical value of T\ =~0. 15. This is consistent 
with the Banfield et al. (1998b) SSI retrievals of T\ =~0.1 and the value of 
t ] (660) = 0.2 found by Sromovsky and Pry (2002) using fits to center-to-limb 
brightness curves from WFPC2 observations. All zones have greater T\ than 
the belts and the vortices, with the exception of the red annulus of Oval BA, 
which was comparable to the zones. The three lowest tj values were found in 
the NEB, NEB1, and GRS center. 

The optical depths of the tropospheric haze (Fig. 20 and Table 6) only deviate 
slightly from the initial value of 7*2 = 4, except for the NEB and NEB1, where 
optical depths increase to 15. Except for these high values in the NEB and 
NEB1, i ~2 fell within the range of 2 > r 2 > 8 reported by West et al. (2004). 
Locations with redder observed colors (e.g. belts, GRS center, Oval 2) had 
slightly higher T 2 overall, while locations with the whitest observed colors 
(e.g. EZ, NTrZ, A5) had slightly lower T 2 . This is consistent with Irwin and 
Dyudina (2002), whose radiative transfer models indicated “that the haze near 
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the tropopause (0.26 bars) is denser above the belt than above the zone.” 


The horizontal variations in our T\ and rjj retrievals are fairly small. When con- 
sidered together, they suggest that Jupiter would be nearly uniformly bright 
in CH4 absorption bands. Observations in these bands, however, show that 
zones and the GRS appear brighter than their surroundings (e.g. Karkoschka, 
1998). This discrepancy between our model retrievals and observations in CH 4 
absorption bands is not surprising, as our data did not include CH 4 filters. 

The optical depths of the tropospheric cloud (Fig. 21 and Table 6) exhibit 
small differences between belts and zones. The NTrZ and NTZ both have 
a greater average r 3 than the belts, but the EZ does not. These results do 
not follow the general consensus that belts contain significantly less opacity 
in the tropospheric cloud deck than zones do. Irwin et al. (2001) suggest, 
however, “that the main source of the near-IR reflectivity variability, observed 
to be anticorrelated with 5-/im brightness, is due to opacity changes of cloud 
layers lying at pressures between 1 and 2 bar, and not to variability of the 
higher altitude ammonia cloud.” Thus, our near-UV and visible data may not 
be sensitive to the altitudes at which the belt/zone opacity differences are 
greatest. Retrieved values for individual locations differ by ~15 within each 
weather region, with most regions in the range of 18 > r 3 > 50. However, 
~50% of the locations in both the NEB and NEB1 have significantly lower 
optical depths (r 3 ~ 10) than all other weather regions. 

An interesting result in the optical depth retrievals is that the tropospheric 
haze is highly positively correlated with the optical depth of the cloud in 
many of the weather regions (Fig. 22 and Table 6), with an overall correlation 
coefficient of 0.147. This correlation is strong in small weather regions (e.g. 
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0.995 for the NEB1 and 0.789 for the GRS center) whose the data cuts have 
inherently limited longitudinal coverage, but it is also equally present in some 
zonal regions as well (e.g. 0.974 for the NEB and 0.818 for the NTBn). Our 
interpretation is that the correlations are indicative of a degeneracy in the 
model, which is easily accounted for by the poor vertical sensitivity of the 
models due to the lack of data at CH 4 absorption bands. Instead of physically 
interpreting the optical depths of the two layers independently, it may be 
necessary to consider differences in their ratios when looking for real horizontal 
variations in the model results. 

4-2 Single Scattering Albedos of Tropospheric Haze Aerosols 

We retrieved the tropospheric haze particle color from a globally representa- 
tive sample of jovian locations (Table 7). Figure 23 shows the mean Wq(\) from 
375-673 nm for all data cuts, with the initial zuo(X) for comparison. Figure 24 
presents the results by weather region. Again, this is the first global analysis 
with multiwavelength retrievals, so we compare our results with the global 
analysis of a single wavelength by Simon-Miller et al. (2001b) and the mul- 
tiwavelength analysis of Oval BA and its immediate surroundings by Perez- 
Hoyos et al. (2009). 

Our values for tuo(410 nm) are in good agreement with the range found by 
Simon-Miller et al. (2001b) of 0.922-1.0 in their analysis of Galileo SSI data. 
Our results are also consistent with zuo values for filters F330W, F435W, and 
F658N reported by Perez-Hoyos et al. (2009) in their study of oval BA (see 
their Fig. 9 and 11). Their ruo(A) results were slightly higher; however, they 
assumed a darker stratospheric haze ( 070=0.95). 
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We subtracted the mean ca 0 (A) (Fig. 23, solid line) and ran a PCA in order 
to explore the spectral shapes that contribute to the variance about the mean 
ro'o(A). We pause to note that all PCA results reported in this work (the 
cc7 0 (A) PCA here and the I/F( A) PCA in §5.3) hold true for two different 
PCA implementations. In the first case, all 1345 data cuts are included in the 
PCA. In the second case, the data cuts in each weather region are averaged 
and the PCA is run on those 14 averages. This demonstrates that our disparate 
number of data cuts between weather regions does not bias the PCA results. 
We present only the results for PCAs of all 1345 data cuts. 

The cuo(A) PCA yielded two significant components, PCI and PC2, which 
account for > 99% of the total variance (Fig. 25, top panel). The higher- 
order PCs (Fig. 25, bottom panel) describe an amount of spectral variance 
much smaller than our confidence in the retrievals due to the degeneracy as 
described in Section 3.4. We also view with suspicion the possibility of PC2 
describing a true physical difference between chromophore populations. PC2 
is primarily an anti-correlation between tu 0 (375) and voq (410-502), and the 
value of tuo(375) was the most susceptible to inaccurate retrievals of the wq 
parameters (§3.4). 

The ■ujq(\) for a location is the sum of the mean zuq(\) and each PC multiplied 
by its amplitude. Figure 26 shows the PCI and PC2 amplitudes by weather 
region. As the amplitude of PCI decreases (i.e. moves in the direction of 
negative amplitudes), it describes increasingly redder particle colors that are 
also darker at all wavelengths. Note that locations with the highest positive 
values of PCI still have a tUo(X) that is red. No location has a PCI value high 
enough to have a blue wo(\) once the mean wo(X) is added back. Positive 
values of PC2 describe particles with a decreased ratio of rc7 0 (375) to w o (410- 
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502) than the mean and vice-versa. The data cuts with the strongest PC2 
contribution to their K7 0 (A) belong to the NEB and NEB1, which is another 
reason to doubt the existence of a physical cause behind PC2. The retrievals of 

cuo(375) for these two regions were affected by degeneracy with Pi and rq 3 

more than most other regions (Figs. 12, 14, 16, and 18). 


5 Discussion 

5.1 Chromophore Particle Colors 

The spectra of the chromophores were quite similar for different areas that 
spanned a large range of observed coloration. All locations had rt7 0 (A) with 
a positive slope, and the slope decreases with increasing wavelength. Our re- 
sults are therefore consistent with the presence of an ubiquitous blue absorber 
in Jupiter’s upper troposphere, which was suggested by Simon-Miller et al. 
(2001a). None of the data cuts had ru 0 (A) with a dramatically different shape 
than the mean. We empirically found a functional form that accurately de- 
scribes the zuo(A) from all retrievals: 

zuo(A) = 1 — A exp(-BA). (7) 

The best-fit parameters for the mean cuo(A) are A = 25.4 and B = 0.0149 
for A in units of nm (Fig. 23, dashed line), with an RMS error of 0.00174. 
Not only does this functional form fit the mean wq(A), but it also fits each 
individual cuo(A) retrieval (fitting parameters A and B separately for each of 
the 1345 spectra) with an average and standard deviation of RMS error of 
0.00195 ± 0.00062. Table 8 contains the average and standard deviation of the 
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parameter fits for all individual wq(X) retrievals by weather region. 


5.2 Chromophore Populations 

For both the Svo(A) and I /Fq(X) PCAs, one chromophore alone can explain 
both the mean spectral shape and PCI by horizontal variations in the chro- 
mophore’s concentration. If there is an underlying physical cause for wo(X) 
PC2, then it may be caused by one or more of three possibilities: (1) spatial 
variations in the vertical aerosol structure that affect //Fo(A), which are not 
explored in the simplified model used here to obtain U7 0 (A); (2) spatial varia- 
tions in the physical properties of the chromophore particles, such as radius, 
shape, temperature, pressure, and UV exposure; or (3) spatial variations in 
the concentration of a second chromophore. 

There are two clues to the interpretation of PCI. The first is the fact that 
locations that appear more red than others in observed color (i.e. more neg- 
ative a ig 0 have zuo(X) that are more red (i.e. more negative - 

and darker at all wavelengths. Although we intuitively expected redder par- 
ticles producing redder observed colors, it was not absolutely necessary. The 
second clue is that the optical depth of the tropospheric cloud is positively 
correlated with the amplitude of PCI with an overall correlation coefficient 
of 0.780 (Fig. 27). This indicates a degeneracy between r 3 and brightness of 
chromophore color and/or a real relationship between the two. If this is merely 
a degeneracy, then why did locations with redder observed color all have red- 
der wo(X) retrievals? On the other hand, if it is a real relationship, why is the 
correlation so strong across all retrievals? 
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There is a physically realistic scenario that could produce a strong correla- 
tion. It is probable that locations with higher 73 would also have increased 
upwelling of NH 3 gas. This might cause an increased mixing ratio of NH 3 ice 
to chromophore particles in the tropospheric haze and perhaps increased coat- 
ing of chromophores with ice in the vicinity of the ammonia cloud, both of 
which would brighten our retrieved cao(A). These possibilities — a degeneracy, 
an increased mixing ratio of NH 3 ice to chromophores, and an increased ice 
coating of chromophores — all lead us to conclude that our zuo(X) results are 
consistent with only one chromophore population. We are unable to distin- 
guish between chromophore populations, because that requires a high level of 
certainty concerning cloud heights, optical depths, and particle radii that are 
too poorly constrained in this current study. 


5.3 Comparison between Statistical Analyses of I/Fq( X) and tuq(X) 

We examined how this 770 (A) PCA compares to an analysis of I /F 0 (X) in 
order to determine the extent to which I /Fq(X) variations can be used as a 
proxy for 770 (A) variations. PCA of the mean-subtracted //Fo (375-673) also 
yielded two significant components with similar amounts of variance to the 
770(A) PCA (Fig. 28, top panel) and higher-order PCs with variances < 1 % 
(Fig. 28, bottom panel). Locations with negative values of PCI have darker 
and redder observed colors than the mean I/Fq(X) (Fig. 4, solid line), and 
those with negative values of PC 2 have particles with an increased ratio of 
375-410-nm absorption to 469-673-nm absorption than the mean. 

Are the largest variations in I/F 0 (X ) correlated with those of 770 (A)? Figure 29 
compares the amplitudes of the I / Fq{ X) PCs and 77q(A) PCs for each jovian 
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location. The strong correlation between the PC amplitudes demonstrates 
the potential usefulness of reflectance spectra for the purpose of determining 
the number of statistically independent chromophores. However, these compo- 
nents may not correspond to identical underlying phenomena, as evidenced by 
the fact that the relationship between the PCs is not one-to-one but contains 
obvious structure. Additionally, the I/Fq{\) PCs have significantly different 
spectral shapes from the zvo(X) PCs. Thus, although the I/Fq{\) PCs do re- 
veal spectral characteristics of a chromophore’s observed color, they do not 
reveal the spectral characteristics of a chromophore’s particle color. 


6 Conclusions 

In this study we derived and statistically analyzed horizontal variations of 
G7o(A) in Jupiter’s troposphere using multispectral HST WFPC2 images. The 
locations in our sample were representative of the typical range of the observed 
jovian color. Our radiative transfer models of I/F yielded the first global 
results for the spectra of chromophore particle colors at multiple near-UV 
and visible wavelengths. An empirical functional form (Eq. 7) fits the G7o(A) 
retrievals very well. 

Although derived values for tx?o(A) from models are dependent upon the as- 
sumed vertical aerosol structure, we do not solve for all parameters of that 
structure. We assumed the pressure levels of the aerosol layers due to the lim- 
ited vertical sensitivity inherent in the data. We derived optical depths and 
chromophore spectral shapes in the near-UV and visible wavelength regime for 
1345 individual planet locations, assuming the chromophores affecting these 
wavelengths were located solely in the tropospheric haze. 
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We characterized the spatial variations in particle color via ti7 0 (A) PCA and 
found that two PCs account for > 99% of the total variance in the spec- 
tral shapes of chromophores. Our results are consistent with the conclusion of 
Siinon-Miller et al. (2001a) that one or two chromophores are necessary to ex- 
plain color variations in the jovian atmosphere. The I / F 0 (X) PCA in this work 
makes a stronger case than Simon-Miller et al. (2001a) for the upper limit to 
the number of chromophores because of our greater spectral coverage. How- 
ever, PC2 in the ru 0 (A) PCA was likely describing spectral shapes that arose 
due to inaccurate retrievals from degeneracy in the radiative transfer model 
parameters — cuo(A) was particularly affected by the choices of the pressures 
and particle radii associated with each aerosol layer. We also find that model- 
ing the physical properties that contribute to the observed I/F( A) provides a 
distinctly different context for characterizing chromophore spectra than PCA 
of I/F 0 { A). We recognize the need for additional multispectral studies that 
sample the continuum below ~600 nm as well as obtain observations in CH4 
absorption bands of various strengths, which will further constrain the radia- 
tive transfer retrievals while making the context more physically appropriate. 

The particle spectra presented in this work represent a large portion of the 
global distribution of chromophore(s) on Jupiter. Yet, we cannot chemically 
identify the chromophore(s) using our derived particle spectra as they contain 
no discrete absorption features — only broad changes across the entire near- 
UV and visible regime. It is possible that an absorption feature may exist in 
the gap in our coverage between 502 nm and 673 nm, which is precisely the 
spectral region identified by Lebofsky and Fegley (1976) where both H2S and 
NH4SH ices develop an absorption feature ~100-nm wide centered at ~600 nm 
when irradiated with UV light (see also Sill, 1973). Moreover, recent work by 
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Sugiyama et al. (2009), Sromovsky and Fry (2010a), and Sromovsky and Fry 
(2010b) indicates NH 4 SH particles may be widespread in Jupiter’s upper tro- 
posphere. Future chromophore studies should include spectral coverage from 
500-700 nm to search for these ice features. We also encourage the continued 
study of chromophore candidates (see West et al., 1986) in the laboratory to 
provide a database of their indices of refraction — as a function of wavelength, 
temperature, pressure, irradiation, etc . — to achieve the ultimate goal of iden- 
tifying the chemical species and the local environmental conditions necessary 
to produce the observed jovian colors that so deeply intrigue us. 
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Table 1. Summary of HST Observations 


UT Date 
(YYYY-MM-DD) 

Spectral 
Cube ID 

WFPC2 

Filter 

Image # 

2008-05-15 

A 

F255W 

ub060105 


A 

F343N 

ub060106 


A 

F375N 

ub060107 


A 

F390N 

ub060108 


A 

F410M 

ub060102 


A 

F410M 

ub06010d 


A 

F437N 

ub060109 


A 

F469N 

ub06010a 


A 

F502N 

ub06010b 


A 

F673N 

ub060101 


A 

F673N 

ub06010c 


B 

F255W 

ub060205 


B 

F343N 

ub060206 


B 

F375N 

ub060207 


B 

F390N 

ub060208 


B 

F4X0M 

ub060202 


B 

F410M 

ub06020d 


B 

F437N 

ub060209 


B 

F469N 

ub06020a 


B 

F502N 

ub06020b 


B 

F673N 

ub060201 


B 

F673N 

ub06020c 

2008-06-28 

C 

F255W 

ub060403 


C 

F343N 

ub060404 


C 

F375N 

ub060405 


C 

F390N 

ub060406 


C 

F410M 

ub060402 


C 

F437N 

ub060407 


C 

F469N 

ub060408 


C 

F502N 

ub060409 


C 

F673N 

ub060401 


C 

F673N 

ub06040c 

2008-07-08 

D 

F255W 

ub060609 


D 

F343N 

ub060604 


D 

F375N 

ub060605 


D 

F390N 

ub060606 


D 

F410M 

ub060602 


D 

F437N 

ub060607 


D 

F469N 

ub060608 


D 

F502N 

ub060603 


D 

F673N 

ub060601 


D 

F673N 

ub06060c 


44 



Table 2. Weather Regions 



Feature 

Planetogr&phic 

System III 

Longitude 


Mo 

Samples per 

Data 

Region 

Type 

Latitude 

Longitude 

Range 

Range 

Range 

Data Cut 

Cuts 

EZ 

zone 

-0.0 ± 2.0 

122.5 ± 19.2 

44.2 db 4.4 

0.21 ± 0.04 

0.27 ± 0.03 

414 ± 68 

12 

NER 

belt 

11.2 ± 0.6 

122.5 ± 19.2 

32.1 ± 7.5 

0.15 ± 0.03 

0.16 ± 0.03 

301 ± 64 

12 

NTYZ 

zone 

21.6 ± 0.4 

122.5 db 19.2 

46.5 ± 6.1 

0.21 ± 0.05 

0.26 db 0.03 

467 ± 69 

12 

NTBs 

belt 

24.5 ± 0.7 

122.5 ± 19.2 

41.4 ± 7.7 

0.18 ± 0.04 

0.22 ± 0.04 

423 db 74 

12 

NTRn 

belt 

27.3 ± 0.8 

122.5 ± 19.2 

45.9 ± 3.2 

0.19 ± 0.03 

0.24 ± 0.03 

477 ± 33 

12 

NTZ 

zone- 

33.9 ± 0.6 

122.5 ± 19.2 

44.6 ± 11.0 

0.20 ± 0.06 

0.21 ± 0.04 

392 ± 78 

12 

GRS bright annulus 

oval 

-22.1 ± 2.1 

111.3 ± 7.6 

7.6 db 3.9 

0.09 ± 0.04 

0.08 ± 0.04 

56 ± 26 

566 

GRS center 

oval 

-23.1 ± 0.8 

111.6 ± 6.4 

4.7 ± 1.1 

0.08 ± 0.03 

0.07 ± 0.04 

52 ± 11 

115 

GRS red annulus 

oval 

-22.2 ± 3.1 

114.4 ± 7.5 

7.7 ± 5.7 

0.11 ± 0.04 

0.11 ± 0.05 

35 ± 17 

180 

BA center 

oval 

-33.4 ± 0.9 

118.4 db 5.0 

2.1 ± 0.7 

0.07 db 0.02 

0.10 ± 0.04 

27 ±7 

98 

BA red annulus 

oval 

-33.4 ± 1.2 

117.9 ± 5.6 

6.0 ± 1.8 

0.08 ± 0.02 

0.10 ± 0.05 

32 ± 9 

190 

Oval 2 

oval 

-24.4 ± 0.6 

129.6 ± 3.2 

1.9 ± 0.4 

0.03 ± 0.01 

0.07 ± 0.03 

23 ± 4 

52 

NEB1 

oval 

15.5 ± 0.5 

89.4 db 4.0 

1.3 ± 0.3 

0.13 ± 0.08 

0.17 ± 0.06 

21 ± 2 

36 

A5 

oval 

-40.7 ±0.7 

92.5 ± 0.0 

2.2 ± 0.5 

0.03 ± 0.00 

0.05 ± 0.00 

30 ± 5 

36 


All values following a “±” sign are standard deviations (not uncertainties) for that region’s data cuts. 



Table 3. HST WFPC2 Filters and Beer’s Law Parameter Fits 


Filter 

Weighted Mean 
Wavelength (nm) 

Absorption 

Parameter 0 

Scattering 

Parameter 0 

F255W 

260.490 

2.7E-07 

1.842 

F343N 

343.416 

9.5E-06 

0.607 

F3T5N 

373.988 

1.1E-05 

0.432 

F390N 

389.319 

1.1E-05 

0.367 

F410M 

408.811 

1.0E-05 

0.301 

F437N 

436.922 

9.5E-05 

0.230 

F469N 

469.439 

2.7E-05 

0.173 

F502N 

501.244 

1.9E-04 

0.133 

F673N 

673.229 

4.3E-03 

0.041 


a The absorption and scattering parameter values are fits to the calculated filter transmissivities (see 
Fig. 6) with this functional form of Beer’s Law; T — exp(—2Ap), where T is (1 - extinction) at the surface 
for a beam of light traveling down to the given pressure and back, p is the pressure, and A is the 
absorption or scattering parameter. 
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Table 4. Variability Flags for Model Parameters 


Parameter 

Stratospheric 

Haze 

Tropospheric 

Haze 

Tropospheric 
Sheet Cloud 

P 

0 

0 

2 

r( 502) 

1 

1 

1 

ro 

0 

0 

0 

too (255) 

1 

3 

0 

too (343) 

1 

3 

0 

too (375) 

0 

1 

0 

too (390) 

0 

1 

0 

too (4 10) 

0 

1 

0 

too (437) 

0 

1 

0 

too (469) 

0 

1 

0 

too (502) 

0 

1 

0 

too (673) 

0 

1 

0 


Flag definitions: 0 denotes fixed parameters, 1 denotes parameters free to vary, 2 denotes a pressure 
parameter fixed to the layer above (creating an infinitesimally thin layer), and 3 denotes zu o parameters 

fixed to the filter redward in the same layer. 
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Table 5. Initial Values for Model Parameters 


Parameter 

Stratospheric 

Haze 

Tropospheric 

Haze 

Tropospheric 
Sheet Cloud 

P [bar] 

0.100 

o 

jfc* 

P 

jo 

© 

0.4, 0.7, 1.0 

r(502) 

0.200 

4.000 

20.00 

ro [fim] 

0.030 

0.900 

2.000 

wo (255) 

0.900 

0.920 

1.000 

O7o (343) 

0.900 

0.920 

1.000 

too (375) 

1.000 

0.920 

1.000 

O7 0 (390) 

1.000 

0.940 

1.000 

070 (4 10) 

1.000 

0.950 

1.000 

o 7o (437) 

1.000 

0.960 

1.000 

O7o(469) 

1.000 

0.975 

1.000 

o?o (502) 

1.000 

0.980 

1.000 

ojq(673) 

1.000 

1.000 

1.000 


The initial model parameters. Units (if any) are indicated in brackets in the parameter column. The initial 
pressures for the tropospheric haze and sheet cloud are both fixed at 0.4 bar, 0.7 bar, and 1.0 bar. Thus, 

there are three initial models. 
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Table 6. Optical Depth Retrievals 


Region 

Stratospheric Haze 
n( 502) 

Tropospheric Haze 
T2(502) 

Tropospheric Cloud 
T 3 (502) 

Correlation Coefficient 
between T 2 and T 3 

All Regions 

0,16 ± 0.04 

4. ± 1. 

30. ± 9. 

0.147 

EZ 

0.22 ± 0.02 

3.5 ± 0.4 

29. ± 6 . 

0.418 

NEB 

0.10 ± 0.06 

10 . db 6 . 

30 ± 20 

0.974 

NTrZ 

0.21 ± 0.04 

3.1 ± 0.6 

35. ± 4. 

0,048 

NTBs 

0.14 ± 0.05 

4.6 ± 0.2 

30. ± 5. 

0.545 

NTBn 

0.17 ± 0.04 

4.2 ± 0.2 

31. ± 6 . 

0.818 

NTZ 

0.20 db 0.02 

4.2 ± 0.3 

44. ± 5. 

-0.839 

GRS bright annulus 

0.17 db 0.03 

4.2 db 0.1 

27. ± 5. 

0.583 

GRS center 

0.11 db 0.04 

4.8 ± 0.4 

28. ± 6 . 

0.789 

GRS red annulus 

0.16 db 0.04 

4.2 ± 0.2 

26. ± 4. 

0.327 

BA center 

0.17 db 0.03 

3.8 ± 0.4 

41. ± 5. 

-0.041 

BA red annulus 

0.20 db 0.02 

4.3 ± 0.2 

40. ± 7. 

0.421 

Oval 2 

0.15 ± 0.03 

4.4 ± 0.2 

26. ± 3. 

-0.005 

NEB1 

0.10 ± 0.04 

9. ± 5. 

30 ± 20 

0.995 

A5 

0.13 zb 0.02 

3.3 ± 0.5 

46. ± 3. 

0.541 
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Table 7. Single Scattering Albedo Retrievals 


Region 

tuo(255) a 

n7o(343) a 

O7o (375) 

o7 0 (390) 

roo(410) 

tz?o(437) 

tuq(469) 

tu 0 (502) 

too (673) 


All Regions 

0.873 

± 

0.025 

0.893 

± 

0.029 

0.902 

± 

0.016 

0.924 

dr 

0.013 

0.946 

dr 

0.010 

0.962 

dr 

0.009 

0.975 

rt 

0.007 

0.984 

dr 

0.005 

0.999 

dr 

0.001 

EZ 

0.865 

± 

0.015 

0.916 

± 

0.009 

0.922 

dr 

0.004 

0.943 

dr 

0.003 

0.960 

dr 

0.002 

0.972 

it 

0.002 

0.980 

dr 

0.002 

0.986 

dr 

0.001 

0.998 

dr 

0.001 

NEB 

0.891 

± 

0.027 

0.878 

± 

0.036 

0.904 

± 

0.029 

0.916 

rt 

0.021 

0.936 

dr 

0.015 

0.947 

dr 

0.010 

0.960 

dr 

0.007 

0.971 

dr 

0.005 

0.996 

± 

0.003 

NTrZ 

0.878 

± 

0.011 

0.933 

dr 

0.011 

0.931 

± 

0.009 

0.948 

dr 

0.006 

0.963 

dr 

0.003 

0.973 

dr 

0.002 

0.981 

dr 

0.001 

0.986 

dr 

0.001 

0.999 

dr 

0.001 

NTBs 

0.863 

db 

0.032 

0.883 

dr 

0.033 

0.888 

dr 

0.016 

0.909 

dr 

0.011 

0.933 

dr 

0.007 

0.950 

± 

0.005 

0.967 

± 

0.004 

0.978 

dr 

0.002 

1.000 

± 

0.000 

NTBn 

0.878 

dr 

0.017 

0.891 

rt 

0.009 

0.897 

dr 

0.008 

0.919 

dr 

0.006 

0.942 

± 

0.004 

0.959 

dr 

0.003 

0.974 

dr 

0.002 

0.983 

dr 

0.002 

1.000 

dr 

0.001 

NTZ 

0.859 

dr 

0.013 

0.907 

± 

0.005 

0.913 

± 

0.008 

0.938 

dr 

0.005 

0.958 

dr 

0.003 

0.972 

dr 

0.002 

0.983 

dr 

0.001 

0.989 

dr 

0.001 

1.000 

dr 

0.000 

GRS bright annulus 

0.863 

dr 

0.012 

0.883 

dr 

0.020 

0.897 

dr 

0.008 

0.922 

dr 

0.006 

0.946 

dr 

0.005 

0.962 

dr 

0.004 

0.975 

dr 

0.003 

0.984 

.dr 

0.002 

0.999 

dr 

0.001 

GRS center 

0.842 

dr 

0.025 

0.860 

dr 

0.028 

0.893 

dr 

0.020 

0.913 

dr 

0.014 

0.936 

dr 

0.009 

0.952 

dr 

0.006 

0.967 

dr 

0.004 

0.978 

dr 

0.003 

0.999 

± 

0.001 

GRS red annulus 

0.878 

± 

0.017 

0.889 

dr 

0.032 

0.899 

dr 

0.011 

0.919 

dr 

0.009 

0.942 

± 

0.006 

0.958 

dr 

0.005 

0.973 

dr 

0.003 

0.982 

dr 

0.002 

0.999 

dr 

0.001 

BA center 

0.874 

dr 

0.010 

0.915 

± 

0.011 

0.925 

dr 

0.008 

0.945 

dr 

0.005 

0.962 

dr 

0.003 

0.974 

± 

0.002 

0.984 

dr 

0.001 

0.990 

dr 

0.001 

0.999 

dr 

0.001 

BA red annulus 

0.903 

dr 

0.013 

0.922 

dr 

0.015 

0.910 

dr 

0.007 

0.931 

rt 

0.005 

0.953 

dr 

0.003 

0.967 

dr 

0.003 

0.979 

dr 

0.001 

0.987 

dr 

0.001 

1.000 

dr 

0.000 

Oval 2 

0.858 

± 

0.025 

0.884 

dr 

0.020 

0.889 

dr 

0.009 

0.913 

dr 

0.006 

0.938 

dr 

0.005 

0.956 

dr 

0.004 

0.971 

dr 

0.002 

0.981 

± 

0.002 

0.999 

± 

0.001 

NEB1 

0.940 

dr 

0.012 

0.924 

zb 

0.034 

0.896 

dr 

0.023 

0.902 

dr 

0.021 

0.922 

rfc 

0.017 

0.936 

dr 

0.010 

0.953 

dr 

0.007 

0.964 

dr 

0.007 

0.995 

dr 

0.002 

A5 

0.893 

dr 

0.006 

0.928 

dr 

0.012 

0.944 

dr 

0.008 

0.958 

± 

0.006 

0.969 

dr 

0.005 

0.978 

dr 

0.004 

0.985 

dr 

0.002 

0.990 

dr 

0.001 

1.000 

dr 

0.000 


° The values reported here for tz7o(255) and wq( 343) are for stratospheric haze particles. 1370 (255) and n7o(343) for tropospheric haze particles are equal to wo (375). 



Table 8. Best-fit Parameters for the Empirical Functional Form for Wq(\) 


Region 

Parameter 

A 

Parameter 

B 

All Regions 

34.1 ± 18.8 

0.0151 ± 0.0020 

EZ 

28.0 ± 8.0 

0.0157 ± 0.0009 

NEB 

4.0 ± 2.9 

0.0093 ± 0.0015 

NTrZ 

17.4 ± 10.4 

0.0145 ± 0.0011 

NTBs 

17.2 ± 8.8 

0.0131 ± 0.0013 

NTBn 

31.7 ± 12.1 

0.0150 ± 0.0014 

NTZ 

79.1 ± 18.9 

0.0182 ± 0.0005 

GRS bright annulus 

39.7 ± 13.4 

0.0157 ± 0.0012 

GRS center 

16.6 ± 12.2 

0.0129 ± 0.0015 

GRS red annulus 

22.6 ± 8.3 

0.0142 ± 0.0011 

BA center 

47.0 ± 22.1 

0.0170 ± 0.0010 

BA red annulus 

44.7 ± 19.5 

0.0164 ± 0.0009 

Oval 2 

31.6 ± 13.1 

0.0149 ± 0.0010 

NEB1 

3.2 ± 1.8 

0.0086 ± 0.0012 

A5 

13.3 ± 2.5 

0.0146 ± 0.0003 


The functional form (Eq. 7) is — X ~ A exp(~~BA) for A in units of nm. 
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Fig. 1. Throughput curves for HST WFPC2 filters (solid black) used in the observa- 
tions, the approximate full-disk albedo spectrum of the data (dotted black) with the 
corrected value for F343N, and the reference full-disk albedo spectrum of Jupiter 
(gray) from Karkoschka (1998). 
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Normalized Total System Response 









Fig. 2. All four mapped data cubes. Cubes A and B are from 15 May 2008, Cube C 


is from 28 June 2008, and Cube D is from 08 July 2008. Locations that are modeled 


in this work are outlined in white. Latitudes are planetographic, and longitudes are 


System III. Colors are approximately true-color. Each color channel is the average 


of the I/F values in the given filters and scaled from 0 to 1: Red=502-673 nm, 


Green=437-469 nm, Blue=255-410 nm. No limb-darkening correction was applied 


which can be seen most readily on the left sides of B and D. These are greenish 


where fio is low in the blue channel images 





110 100 90 80 70 

System lit Longitude 


Fig. 3. An example data cut (areas outlined in white) selected from the EZ in Cube 
A is shown with the same color channels as in Fig. 2. The samples in the data cut 
that were selected for modeling are the averages of pixels within 1° bins in longitude, 
which are the size of the minor tick marks in this figure. 



Fig. 4. The mean and standard deviation of the I /Fq spectra for all 1345 data cuts. 
The reference full-disk albedo spectrum from Karkoschka (1998) lies below the I/Fo 
spectra because it does not account for limb darkening. 


54 




Wavelength (nm) 


Fig. 5. The mean (solid black) and standard deviation (gray) of I/Fq(X) for all data 
cuts separated by weather region. The dashed black line is the mean for all data cuts 
from all weather regions. All panels have the same abscissa and ordinate ranges. 






















Methane Gas Absorption Rayleigh Scattering 




Fig. 6. Two-way transmissivities for HST WFPC2 filters (solid curves) assuming 
a nadir view and no aerosol opacity. The left panel contains gas absorption trans- 
missivities calculated from the Karkoschka (1998) methane absorption coefficients. 
Curves for all but the F673N filter are overlapping at a transmissivity value of unity 
from space down to ^10 bars. The right panel contains two-way transmissivities for 
Rayleigh scattering. The curves progress from F255W to F673N in order of wave- 
length. The dotted curves in both plots are the Beer’s Law fits used in the radiative 
transfer code (see Table 3). Most fits coincide with the actual transmissivities at 
this resolution. 


56 



Pressure Increasing 



Stratospheric Haze 
Optical Depth tj 

(Colored in ultraviolet. White in visible) 




Tropospheric Sheet Cloud 
Optical Depth X 3 


(White in all wavelengths) 


0.1 bar Fixed Boundary in All Models 


-p 0.4 bar 
-- 0.7 bar 
— 1.0 bar 


Three Models with 
Fixed Pressures 


Fig. 7. Schematic of the vertical structure models used in this work. 
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Data and Best-Fit Model (^=0.62 ) 


255 nm: 1 Image, 48 Samples 343 nm: 1 Image, 48 Samples 375 nm: 1 1mage, 48 Samples 



0,86 0.88 0,90 0.92 0.94 0.86 0.98 1.00 0.80 0.92 0.94 0.96 0.98 1.00 0.90 0.82 0.94 0.96 0.98 1.00 

390 nm: 1 image, 48 Samples 410 nm: 2 Images, 96 Samples 437 nm: 1 Image, 48 Samples 



0.88 0.90 0.92 0.94 0.98 0.98 1,00 0.80 0.86 0.90 0.95 1.00 0.85 0.90 0.95 1.00 



Model Parameters 


w" 0.1 

l Layer 1 

P,=0.100 bar 

r r =0 32 

r t =0.03//m % 

I 

Cl 

t Layer 2 
•• — 

P J? =0.700 bar 

3.52 

r 2 =0.90jarri ;! 

1.0 

EL. , _ 

Layer 3 

P 3 =0.700 bar 

r 3 = 28.46 

r 3 =2.00^m J 



Fig. 8. Example data cut and best-fit model using the ^-minimization routine (see 
Fig. 9). The data cut (Fig. 3, outlined areas) is from the EZ in Cube A. The samples 
in the data cut that were selected for modeling (black points in top plots) are the 
averages of pixels within 1° bins in longitude. The wavelength, and the number of 
images and samples at that wavelength, are noted above each plot of I/F versus fi . 
The error bars are the estimated errors used to calculate y 2 . The green points are 
the radiative transfer results generated by the model parameters at the bottom of 
the figure. 
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0.974 0.978 0.982 0.986 0.9791 0.9827 0.9864 0.9901 0.988 0.995 0.997 1.001 


Parameter Value 


Fig. 9. Example convergence of parameters using the ^-minimization routine. The 
data cut (Fig. 3, outlined areas) is from the EZ in Cube A. The diamonds are from 
models saved during Phases A and B, and the asterisks are the best-fit model at 
the completion of Phase C (see Fig. 8). The dotted lines mark the initial parameter 
value (Table 5). 
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0.970 0.915 0.956 0.957 0.978 


0.885 


X 3 



10 20 30 40 50 


G3 0 (375) 



0.857 0.892 0.928 0.963 


G3 0 (437) 



0.931 0.949 0.967 0.985 



Parameter Value in Synthetic Data 


Fig. 10. Retrievals from synthetic data with added Gaussian noise (a — 0.005 [I/F]). 
The difference between the parameter values used to generate the synthetic data 
and those retrieved are marked at values of 0% (solid line), 1% (dotted line), and 
10% (dashed line). 
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T3| 

<L> 

> 



0.7S 0.S0 0.82 0.54 0.55 0.88 0.900.92 0.80 0.85 0.90 0.95 0.861 0.887 0.913 0.939 



0.880 0.906 0.931 0.957 0.91 0.93 0.95 0.97 0.926 0.943 0.961 0.976 



Parameter Value in Synthetic Data 


Fig. 11. Retrievals from synthetic data with the stratospheric/tropospheric haze 
interface (Pi) pressure set to 0.2 bar. The retrievals assumed Pi =0.1 bar. The 
difference between the true parameter values and those retrieved are marked at 
values of 0% (solid line), 1% (dotted line), and 10% (dashed line). 
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EZ NEB NTrZ 



400 500 600 700 400 500 600 700 400 500 600 700 


NTBs NTBn NTZ 



GRS bright onnulus GRS center GRS red onnulus 



BA center BA red onnulus Oval 2 



400 500 600 700 400 500 600 700 400 500 600 700 


NEB1 A5 



400 500 600 700 400 500 600 700 

Wavelength (nm) 

Fig. 12. Retrievals from synthetic data with the stratospheric/ tropospheric haze 
interface (Pi) pressure set to 0.2 bar. The retrievals assumed Pi— 0.1 bar. The solid 
lines are the ruo(A) values of the synthetic data, and the gray lines are the retrievals. 
There are three models (P3 = [0.4, 0.7, 1.0] bar), and all three were used for these 
retrievals. In most of these plots the three gray lines overlap, because all three 
models usually converged to similar zuq(A) curves. 
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Parameter Value in Synthetic Data 


Fig. 13. The same as Fig. 11, except for retrievals from synthetic data with the 
radius of stratospheric haze particles set to 0.1 fi m. The retrievals assumed 0.03 ^m 
particles. 
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Fig. 14, The same as Fig. 12, except for retrievals from synthetic data with the 
radius of stratospheric haze particles set to 0.1 /xm. The retrievals assumed 0.03 fim 
particles. 
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®0 (343) 
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03 0 (673) 




Parameter Value in Synthetic Data 


Fig. 15. The same as Fig. 11, except for retrievals from synthetic data with the 
radius of tropospheric haze particles set to 2.0 fi m. The retrievals assumed 0.9 (x m 
particles. 
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BA center BA red annulus Oval 2 
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Fig. 16. The same as Fig. 12, except for retrievals from synthetic data with the 
radius of tropospheric haze particles set to 2.0 /mi. The retrievals assumed 0.9 /un 
particles. 
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Parameter Value in Synthetic Data 


Fig. 17. The same as Fig. 11, except for retrievals from synthetic data with the 
radius of tropospheric cloud particles set to 4.0 fim. The retrievals assumed 2.0 fim 
particles. 
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Fig. 18. The same as Fig. 12, except for retrievals from synthetic data with the 
radius of tropospheric cloud particles set to 4.0 fim. The retrievals assumed 2.0 fim 
particles. 
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NTrZ NTBs NTBn 





Fig. 19. The optical depth at 502 nm of the stratospheric haze (ri). All panels have 
the same abscissa and ordinate ranges. The dashed line is located at value of r\ of 
the initial model provided to the ^-minimization routine. The number AT denotes 
the number of data cuts modeled in that weather region. The bin size is 0.01 optical 
depths. 
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Fig. 20. The same as Fig. 19, except for the optical depth of the tropospheric haze 
( 72 ). The bin size is 0.5 optical depths. 
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Fig. 21. The same as Fig. 19, except for the optical depth of the tropospheric cloud 
( 73 ). The bin size is 1.5 optical depths. 
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Fig. 23. The mean (solid black) and standard deviation (gray) of the retrieved 
single scattering albedo spectra (z*7o(A)) for all data cuts. The dashed black line is 
a functional fit to the mean with the form given in Eq. 7. The dotted black line is 
the initial model provided to the x 2 minimization routine. 
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Fig. 24. The mean (solid black) and standard deviation (gray) of the retrieved 070(A) 
for all data cuts separated by weather region. The dashed black line is the mean 
for all data cuts from all weather regions and is repeated in each panel for ease of 
comparison between regions. All panels have the same abscissa and ordinate ranges. 
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Fig. 25. Principal components (PCs) from a principal component analysis (PC A) 
of the derived single scattering albedo spectra for all data cuts. The number to the 
right of each PC denotes the percent of the total variance in the spectra described 
by that PC. The top panel shows PCI and PC2. The bottom panel contains PCs 
3-7 with a different ordinate scale. The gray dotted lines are at the same coefficient 
values in both panels for comparison between the two ordinate scales. The zero 
value is indicated with a black dotted line. For each PC, coefficients with the same 
sign are positively correlated, and coefficients with opposite signs are negatively 
correlated. 
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53 - 0,002 0,000 0,002 - 0.002 0.000 0.002 - 0.002 0.000 0.002 

O 
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-0,002 0.000 0.002 -0.002 0.000 0.002 -0.002 0.000 0X02 


Oval 2 NEB1 A5 



-0X02 0,000 0.002 -0.002 0.000 0.002 -0.002 0.000 0.002 


PC Amplitude (arbitrary) 

Fig. 26. Amplitudes of PCI (black) and PC2 (gray) for each weather region. All 
panels have the same abscissa and ordinate ranges. The dashed line is located at 
an amplitude of zero. The number N denotes the number of data cuts modeled in 
that weather region. The bin size is 0.0001. 
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Fig. 27 . The optical depth of the tropospheric cloud (73) versus the amplitude of 
PCI from the u 70(A) PCA. All panels have the same abscissa and ordinate ranges. 
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Fig. 29. Amplitudes of PCs from the I /Fq(\) PCA are compared with those from 
the o7o(A) PCA. The top panel shows the amplitudes of PCI, and the bottom panel 
shows the amplitudes of PC2. 




